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Buckling  of  Unidirectional  Graphite/Epoxy  Composite 
Plates  at  Low  Temperatures 

PIYUSH  K.  DUTTA,  DAVID  HUI,  AND  YVONNE  M.  TRAYNHAM 


INTRODUCTION 

Buckling  and  postbuckling  behavior  of  composite 
material  plates  has  been  thoroughly  investigated  because 
of  their  many  aerospace,  automotive,  marine  and  sporting 
goods  applications  (Singer  1985,  Leissa  1985,  Chen 
and  Chen  1987,  Biswas  1976,  Chen  et  al.  1985,  Flaggs 
and  Vinson  1978,  Jones  1975).  Surprisingly,  such  be¬ 
havior  under  cold  environments  has  received  only  scant 
attention.  Even  in  the  simple  case  of  unidirectionally 
graphite-epoxy  material  plates,  no  buckling  studies 
have  been  reported  in  the  literature. 

The  objective  of  this  investigation  is  to  study  the 
effect  of  low  temperatures  on  the  unidirectional  graph¬ 
ite  epoxy  laminate  plates  in  the  buckling  mode  when 
they  are  constrained  and  loaded  in  the  direction  of  the 
fibers.  The  current  work  is  particularly  interesting  since 
such  graphite/epoxy  plates,  due  to  their  negative  ther¬ 
mal  expansion  coefficients,  expand  in  the  fiber  direc¬ 
tion  under  cooling,  whereas  in  most  conventional  mate¬ 
rials,  lowering  the  temperature  will  cause  contraction. 
Under  in-plane  constraint,  such  graphite-epoxy  plates 
may  develop  out-of-plane  bending  under  cooling  even 
without  any  mechanical  load.  Thus,  upon  bending,  the 
in-plane  load  canying  capacity  of  the  plate  may  be 
significantly  reduced.  The  same  observation  can  be 
made  on  thermal  heating  of  other  conventional  materi¬ 
als  having  positive  thermal  expansion  coefficients. 
Normally ,  a  plate  made  of  conventional  material  becomes 
stiffer  under  thermal  cooling  and  the  critical  buckling 
load  is  expected  to  be  higher.  But,  in  the  case  of 
unidirectional  graphite-epoxy  plates  undercooling,  there 
are  two  opposing  effects:  the  structure  becomes  stiffer 
(advantage)  and,  in  the  presence  of  in-plane  constraint, 
the  structure  will  bend  due  to  expansion  (disadvantage). 
Thus,  it  is  of  interest  to  determine  which  effect  will 
dominate  in  the  buckling  problems.  The  analysis  pre¬ 
sented  here  is  novel  in  that  the  plate  bending  due  to 
thermal  cooling  is  treated  as  an  initial  geometric  imper¬ 
fection  in  the  buckling  and  postbuckling  investigation. 


The  basic  findings  of  the  unidirectionally  laminated 
plate  will  enhance  the  understanding  of  the  buckling 
behavior  of  composite  plates  under  thermal  environ¬ 
ments. 

In  the  current  study,  a  number  of  unidirectional 
graphite-epoxy  laminated  plates  were  tested  at  room 
and  low  temperature  under  uniaxial  compressive  load. 
To  determine  the  bending  strains,  the  plates  were  instru¬ 
mented  with  strain  gauges  mounted  on  both  sides  of 
their  surfaces  at  the  center.  The  load  was  applied  in  a 
servohydraulic  testing  machine,  and  strains  were  re¬ 
corded  using  a  data  acquisition  system. 

BUCKLING  OF  PLATES  UNDER 
THERMO-MECHANICAL  LOADS 

The  bending  of  the  plate  due  to  thermal  load  is 
treated  as  an  initial  geometric  imperfection  having  the 
same  shape  as  the  buckling  mode.  Such  imperfection 
amplitude  can  be  obtained  as  a  function  of  the  tempera¬ 
ture.  The  postbuckling  behavior  is  modeled  using 
Koiter’s  theory  of  elastic  stability.  In  the  analysis,  the 
plate  at  room  temperature  is  assumed  to  be  perfectly  flat 
(zero  imperfection  amplitude)  and,  prior  to  the  applica¬ 
tion  of  mechanical  loads,  the  imperfection  is  assumed  to 
be  due  to  thermal  loading  only. 

The  thermal  expansion  coefficient  of  unidirectional 
graphite-epoxy  plates  in  the  fiber  direction  is  known  to 
be  negative  so  that  that  the  material  expands  upon 
cooling.  Tsai  and  Hahn  (1980)  reported  that  a  typical 
value  for  AS/3501  is-0.3x  10'^per°C.Thevaluesvary 
from  zero  to  -3.0  x  10  per  °C,  and  typical  values  were 
reported  by  Freeman  and  Campbell  (1972),  Lord  and 
Dutta  (1988)  and  Rondeau  et  al.  ( 1 988).  The  experimen¬ 
tal  value  of  the  thermal  expansion  coefficient  of  the 
current  set  of  graphite/epoxy  composites  determined  by 
a  Perkin-Elmer  Thermo  Mechanical  Analysis  System 
was  -0.44  X  1 0“^  per  °C  (Dutta  1991). 

The  rectangular  plates  were  simply  supported  at  the 


(4) 


two  loaded  opposite  edges  and  free  in  the  remaining  two  6  =  5  V(  1  -  A.) 

unloaded  edges.  Such  plates  are  known  as  wide  beams. 

It  should  be  noted  that  the  buckling  mode  of  a  wide  It  should  be  noted  that  the  postbuckling  theory  is  based 
beam  involves  waviness  both  in  the  loading  jr-direction  on  an  asymptotic  analysis,  which  is  valid  only  for  a 

and  in  the  in-plane  y-direction  (Fig.  1 ),  whereas  the  sufficiently  small  amplitude  of  the  imperfection  (say  8/ 

buckling  mode  of  an  Euler  column  involves  waviness  h  <  1 .5). 

only  in  the  j:-direction.  Donnell  (1976)  found  that  the  The  thermal  strain  &j.  in  the  jc-direction  of  the  plate  is 

buckling  load  using  the  wide  beam  (plate)  theory  is  only  given  by 

a  few  percent  higher  than  that  using  the  Euler  buckling 

formula.  In  the  present  analysis,  it  is  assumed  that  the  £.p  =  AZ./L  =  (a)  (AT)  (5) 

postbuckling  behavior  of  wide  beam  can  also  be  predicted 

using  the  Euler  column  theory.  Thus,  in  our  analysis,  the  where  a  is  the  coefficient  of  thermal  expansion  in  the 

influence  of  waviness  in  the  y-direction  is  neglected  fiber  direction  (negative  means  that  the  plate  will  ex¬ 
it  is  well  known  that  the  axial  buckling  load  of  a  pand  upon  cooling)  and  AT  is  the  current  tempierature 

simply  supported  Euler  column  is  minus  the  reference  temperature  (room  temperature). 

The  thermal  buckling  problem  occurs  only  when  (a)  (AT) 
=  tO-  EIIL^  ( 1 )  is  positive  with  in-plane  constraint. 

If  the  deflection  profile  of  the  simply  supported  plate 
where  E  is  elastic  modulus  in  the  jr-(fiber)  direction,  /,  due  to  thermal  loading  alone  is  assumed  to  be  of  the 

the  moment  of  inertia  (=  Bh~ll2),  L,  the  plate  length,  B,  sinusoidal  shape, 
the  width,  and  h,  the  plate  thickness.  _ 

The  postbuckling  behavior  of  an  Euler  column  was  w  =  8  sin  (it  x/L)  (6) 

analyzed  by  Koiter  ( 1 945)  and  Hui  and  Hansen  ( 1 980), 

among  others.  Using  the  terminology  employed  by  where  w  is  the  out-of-plane  deflection,  x  is  measured 

Budiansky  and  Hutchinson  (1964)  and  Hui  (1984),  the  from  the  loaded  edge  and  8  is  the  amplitude  of  the  “im- 

postbuckling  behavior  is  described  by  perfection”  due  to  thermal  loading.  The  arc  length  L  -t- 

AL  can  be  obtained  from  (Fig.  2): 

{bWhf  +  (1  -  [  F/Pcr])  (8//I)  =  (P/Fcr)  (8//1)  (2)  ^ 

where  b  is  the  postbuckling  coefficient,  8  is  the  ampli-  L  +  AT-  -  I  [l  +  (dw/riv)]  dx- 

tude  of  the  initial  geometric  imperfection,  8  is  the  am-  jx  =  0 

plitude  of  the  out-of-plane  deflection  (measured  from  Substituting  the  deflection  w  into  the  arc  length  expres- 

the  imperfect  position)  and  F  is  the  applied  compressive  sion  and  imegrating  and  solving  for  the  imperfection 

load.  From  Hui  and  Hansen  ( 1 980),  it  can  be  shown  that  amplitude  8,  one  obtains 

6  =  (ji2/8)  (/j/T-)2  (3)  8/;i  =  [2/(jt/i)]  (Z-)  (AZ./Z-)'^  (8) 

which  is  essentially  zero  for  /i  <  <  T-  in  plates.  Defining  where  h  is  the  thickness  of  the  laminated  plate.  From  the 

X  =  PIP^^,  it  follows  that  thermal  strain  (eq  5),  the  imperfection  amplitude  can 


Figure  1 .  Reference  coordinates  of  the  buckling  problem. 
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Figure  2.  Deflection  profile  of  the  plate  from  thermal 
stress. 

also  be  written  as 

l/h  =  [2/(rt/i)]  (L)[a(A7)]  .  (9) 


From  these  strains,  the  deflection  can  be  obtained 
using  thecurvature  expression  (Roark  and  Young  1975): 

IIR  =  MI(EI)  (12) 

Since,  in  bending  of  a  beam  M  =  (El){(Pwldx^),  it  fol¬ 
lows  that 

llR  =  ePw/dx2  (13) 

where  R  is  the  radius  of  curvature  of  the  beam  and  M  is 
the  bending  moment.  Note  that  both  R  and  M  are 
functions  of  x,  and  the  product  R,  and  A/  is  a  constant.  If 
the  deflection  shape  is  assumed  to  be  a  half  sine  wave, 

w  =  (8 -t- 5)  sin  (Jur/L)  .  (14) 


Solving  for  the  load  ratio  X,  and  using  eq  9,  one  obtains  Substituting  the  deflection  into  eq  12  and  13,  one  obtains 


X= - ^ -  .  S  +  5  =  L2/(/?n2)  .  (15) 

(5//i)  +  [2L/(jt/i)]  [a(A7-)]'^  (10) 

From  Figure  3,  the  angle  6  is  defined  to  be 

In  order  to  detect  bending  of  the  plate,  the  strains  are 
measuredatthecenteroftheplateonbothsurfaces(Fig. 

3).  The  measured  strains  at  aa  and  66  are,  respectively,  ~/f  -(6/2)  ~^  +  (6/2) 


~  -  1  <  0  compression  so  that 

(11) 

Er  =  (V-^c^  - 1  >  0  tension  R  =  (6/2)  lis^/s^)  +  (55/5^)]  /  [(•^b/'^c^  ~ 


where  s^,  Sg  and  s^^  denote  the  arc  length  aa,  66  and  cc. 

In  the  case  of  th^ermal  loading  only  (no  mechanical 
load),  the  measured  strains  are  denoted  by  e^j^nd  63 • 
Since  the  tensile  and  compressive  modulus  are 
different,  will  not  be  the  same  as  Er  .  Likewise,  in 
general,  is  not  the  same  as  Eg. 


Figure  3.  Bending  of  the  plate. 


(17) 

Substituting  sjs^.  and  s^/s^  from  eq  1 1  into  eq  1 5  and 
17,  one  obtains 

(8  +  8)/6  =  ^  ~  .  (18) 
2  +  £a  +  £b 

Since  £^  +  £g  <  <  2,  no  significant  error  will  occur  from 
their  omission  in  the  denominator.  The  deflection  8/6 
can  then  be  obtained  from  the  offsetted  strains  £g  -  Eg 
and  £^  -  E^  such  that 

8/6  =  [L/(6ji)]2  [(Eg  -  £0  )  -  (E^  -  gO  )]  (19) 

where  and  Eb  are  the  strains  due  to  thermal  load¬ 
ing  only. 


EXPERIMENTAL  METHOD 

The  experimental  buckling  behav  ior  of  the  unidirec¬ 
tional  graphite  fiber/epoxy  plates  was  examined  using 
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Table  1.  Material  specifications  of  the  test  com¬ 
posite  panels. 


Manufacturer 

Fiberite  Corporation 

Fiber  type 

T-300  graphite 

Resin 

Fiberite  974  epoxy 

No.  of  plys/lay  up 

Eight/unidirectional 

Thickness/ply 

0.136  mm  (0.00537  in.) 

Total  thickness 

0.965  mm  (0.038  in.) 

Density 

1522.4  kg/m3  (0.055  lb/in.3) 

Fiber  vol  % 

55 

Major  Poisson’s  ratio 

0.25 

Longitudinal  compr.  modulus 

147.6GPa(21.4x  10<i|bf/in.2) 

Longitudinal  tension  modulus 

169.7  GPa  (24.6  x  lO®  lbf/in.2) 

Transverse  tension  modulus 

1 0. 1 4  G  Pa  ( 1 .47  X  1 0*  lbf/in.2) 

Shear  modulus 

4.62  GPa  (0.67  x  10*  Ibf/in.^) 

Longitudinal  tensile  strength 

2.621  GPa  (380  X  103  lbf/in.2) 

Longitudinal  compr.  strength 

1 .034  GPa  ( 1 50  X  1 03  lbf/in.2) 

Transverse  tensile  strength 

60.0  MPa  (8.7  x  103  lbf/in.2) 

a  servohydraulic  testing  machine  to  apply  a  quasi-static 
load  in  the  fiberdirection  (0°).  The  frame  of  the  machine 
is  enclosed  in  an  environmental  chamber  connected  to 
a  refrigeration  unit  to  test  at  low  temperatures.  The 
plates  were  supported  at  the  edges  using  bearing- 
mounted  components  to  allow  for  free  rotation  at  the 
edges.  At  the  center  of  each  plate,  rosette  strain  gauges 
were  mounted  “back-to-back”  on  either  side.  The  out¬ 
puts  of  the  six  strain  gauges  were  passed  through  a 
Wheatstone  bridge,  while  the  applied  load  and  displace¬ 
ment  were  monitored  and  recorded  using  a  datalogger 
system  and  then  sent  directly  to  a  personal  computer. 

The  composite  plates  tested  in  the  experiment  were 
manufactured  by  Fiberite  Corporation  as  eight-ply  uni¬ 
directional  T-300  graphite  fibers  in  a  974  epoxy  resin 
matrix.  The  material  specifications  are  listed  in  Table  1 . 

The  plates  were  carefully  cut  to  avoid  heating  the 
edges  to  127  x  254  mm  (5  x  10  in.)  using  a  diamond 
blade  on  a  band  saw.  The  plates  were  then  examined 
with  a  dial  gauge  to  determine  the  flatness  and  initial 


Load 


imperfections  of  each  sample.  These  readings  indicated 
that  the  plates  were  flat  with  variations  in  the  thickness 
of  up  to  0.125  mm  (0.005  in.).  Rosette  strain  gauges 
were  then  mounted  on  the  surface  of  each  plate  on  both 
sides  in  a  back-to-back  configuration. 

The  experimental  setup  is  illustrated  in  Figure  4.  The 
system  consisted  of  the  plate  mounted  on  the  bearings 
which  were  supported  by  specially  designed  fixtures  on 
the  top  and  the  bottom,  and,  for  negligible  support,  by 
foam-backed  plexiglass  plates  on  the  sides  that  were 
made  to  fit  against  the  sides  of  the  environmental  chamber. 


Figure  5.  Bearing-mounted  simple  support. 
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Figure  6.  Photograph  of  the  hearing-mounted  simple  support. 


The  design  of  the  bearing  mounted  supports  is  shown 
in  Figures  5  and  6.  A  19-mm  (3/4-in.)  Teflon  rod 
machined  to  create  a  1.3-mm  (0.05-in.)  slot  and  a  1.6- 
mm-  (0.063-in.-)  diam.  hole  along  the  entire  length  and 
then  cut  into  25.4-mm  ( 1-in.)  cylinders  was  used  as  the 
body  of  the  support.  A  semicircular  identation  was 
made  to  a  depth  of  3.2  mm  (0. 125  in.)  to  accommodate 
the  12.7-mm  (0.5-in.)  stainless  steel  ball  bearing.  The  bear¬ 
ings  were  attached  to  the  Teflon  using  a  thin  coat  of  adhesive. 

The  loading  frames  on  the  top  and  bottom  were  made 


from  19-mm  (3/4-in.)  aluminum  plates  with  a  25.4-mm 
( 1-in.)  diameter  hole  in  the  center  tapped  to  14  turns  per 
inch  to  fit  onto  the  testing  machine  thread.  The  plates 
served  as  a  flat  loading  surface  through  the  ball  bearings 
during  the  test.  The  load  was  applied  at  the  rate  of 0.625- 
mm  (0.025-in.)  per  minute  cross  head  speed. 

The  servohydraulic  testing  machine  with  its  envi¬ 
ronmental  chamber  attached  to  a  refrigeration  unit  is 
shown  in  Figure  7.  Insulation  was  placed  at  the  openings 
around  the  upper  and  lower  pistons  of  the  testing  ma- 


Figure  7.  MTS  servohydraulk  test  machine  w  ith  environment  chamber. 


Figure  8.  Schematic  of  test  and  data  collection  system. 


chine  and  also  at  a  side  port  for  the  strain  gauge  wires. 
This  allowed  for  the  strain  gauge  instrumentation  and 
the  load  cell  to  be  protected  at  ambient  temperature. 

The  test  procedure  required  that  the  specimen  be 
perfectly  aligned.  This  was  achieved  by  first  attempting 
to  visually  set  up  the  specimen  supports  along  the 
guidelines  of  the  fixtures,  adjusting  the  strain  gauge 
electrical  instrumentation,  applying  a  small  pre-load, 
and  then  carefully  adjusting  the  specimen  by  the  strain 
gauge  output  monitored  by  the  computer. 

After  alignment  of  the  specimen,  the  load  was  ap¬ 
plied  until  the  specimen  displayed  a  significant  out-of¬ 
plane  displacement.  The  load  was  then  removed  from 
the  plate  and  the  data  were  downloaded  to  a  computer 
for  analysis.  A  schematic  of  the  instrumentation  and 
data  collection  system  is  shown  in  Figure  8. 


EXPERIMENTAL  RESULTS 

The  experimental  results  indicate  some  interesting 
behavior  in  the  buckling  and  postbuckling  modes  for  the 
composite  plates.  The  effect  of  lower  temperature  on 
this  behavior  is  apparent  upon  comparison  of  the  results 
to  those  from  tests  done  at  ambient  temperature. 

Table  2  gives  the  summary  of  the  results.  The  data 
from  a  typical  test  are  shown  in  Figure  9  in  which  the 
applied  stress  and  resultant  axial  strains  are  plotted  with 
respect  to  time.  It  can  be  seen  that  the  stress  increases 
proportionally  with  the  strain  until  the  point  of  bifurca- 


Table  2.  Experimental  buckling  results. 


Test 

Plate 

Temperature 

(X) 

Maximum 

stress 

MPatll^inTt 

Buckling 

stress 

MPa  (tbflin.^) 

C2AR 

C2 

21.0 

1.692(245.38) 

1.645  (238.54) 

R4DR 

R4 

21.0 

1.664(241.33) 

1.418(205.64) 

C4AR 

C4 

21.0 

1.343(194.67) 

0.529  (76.77) 

C3AR 

C3 

21.0 

2.231  (323.55) 

1.730(250.92) 

R3AC 

R3 

-6.1 

1.210(175.52) 

0.416(60.3) 

C5BC 

C5 

-9.6 

1.314  (I90..56) 

1.229(178.26) 

C5DC 

C5 

-11.2 

1.116(161.75) 

0.492  (71.28) 

R4BC 

R4 

-11.3 

1.645(238.58) 

1.277(185.10) 

R4AC 

R4 

-13.3 

1.626(235.79) 

1.182(171.37) 

R2CC 

R2 

-14.0 

1.370(198.77) 

0.567  (82.25) 

R3BC 

R3 
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Figure  9.  Stress  and  strain  history  of  buckling  load. 


tion  of  the  outputs  from  the  two  back-to-back  strain 
gauges,  A  and  B,  on  opposite  surfaces  at  the  center  of 
the  plate.  This  point  of  bifurcation  indicates  that  the 
specimen  has  buckled,  and  the  out-of-plane  displace¬ 
ment  occurred  such  that  one  side  of  the  plate  had  gone 
in  compression  as  the  other  side  went  in  tension.  How¬ 
ever,  the  plate  continues  to  support  additional  loading  in 
the  post-buckling  range  as  the  strain  continues  to  increase. 
The  highest  load  sustained  and  the  buckling  load  found 
from  the  point  of  bifurcation  of  the  axial  strain  gauges 
are  listed  in  Table  2. 

The  axial  strain  at  four  different  temperatures  given 
by  gauges  A  and  B  vs  stress  is  shown  graphically  in 
Figure  10.  A  review  of  these  graphs  demonstrates  a 
marked  difference  in  the  buckling  and  postbuckling 
paths  as  the  temperature  is  decreased.  Figure  10a  shows 
that,  following  a  linear  elastic  deformation,  at  room  tem¬ 
perature  the  plate  buckles  at  a  well-defined  critical  load 
^cr  (approximately  1 .2  MPa  (175  Ibf/in.^).  But  at  lower 
temperatures  the  sharp,  well-defined  buckling  load  is 
absent:  the  curves  of  buckling  path  are  progressively 
flatter  from  -10°  to  -19°C  as  shown  in  Figures  lOb 


through  lOd.  A  series  of  smooth  curves  drawn  through 
the  data  points  of  each  of  the  above  figures  are  superim¬ 
posed  and  shown  in  Figure  lOe.  It  is  clear  from  Figure 
1  Oe  that  lower  temperature  represents  a  higher  degree  of 
initial  geometric  imperfection,  as  discussed  before. 

ITie  stress  vs  axial  strain  is  plotted  for  small  loads 
(before  buckling)  in  Figure  1 1 ,  for  a  representative  plate 
using  a  regression  analysis  to  confirm  that  the  modulus 
of  elasticity  is  consistent  with  the  specifications  pro¬ 
vided  by  the  manufacturer.  These  graphs  serve  only  as 
an  indicator,  as  these  plates  do  not  have  a  suitable 
geometry  to  determine  the  compressive  moduli.  How¬ 
ever,  it  does  show  the  characteristic  stress-strain  behav¬ 
ior.  which  was  within  the  range  of  specifications  of  the 
laminate  plate. 

In  Table  2  the  data  are  arranged  in  order  of  ambient 
to  decreasing  temperatures  to  see  the  effect  on  the  load 
capacity.  Four  basic  ranges  of  temperature,  21°,  -10°, 
-15°,  and  -20°C,  are  grouped.  The  buckling  stresses 
were  determined  by  an  estimate  of  the  stress  at  the  point 
of  departure  from  linearity.  The  maximum  stress  is 
usually  the  stress  at  which  the  test  was  terminated. 
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Figure  10.  Buckling  and  postbuckling  paths  at  different  temperatures. 
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Figure  11.  Axial  stress-strain  data  of  a  typical  panel 
before  buckling. 


DISCUSSION  OF  RESULTS 

The  focus  of  this  study  is  to  determine  experimen¬ 
tally  the  effect  of  low  temperature  on  the  buckling 
behavior  of  graphite  fiber-epoxy  composite.  The  unidi¬ 
rectional  simply-supported  plate  was  chosen,  so  that 
this  effect  could  be  isolated  from  other  interlaminar 
stresses  for  examination.  The  results  show  interesting 
characteristics  that  are  discussed  here. 

In  testing  the  specimen  the  bearing  mounts  described 
above  sucessfully  provided  simple  support  on  the  top 
and  bottom  and  the  boundary  conditions  on  the  sides 
provided  free  edges.  The  plates  did  buckle,  having 
simple  support  in  the  axial  direction  as  an  Euler  column, 
as  described  in  the  theory  section. 

The  results  show  very  significant  changes  in  the 
buckling  and  postbuckling  behavior  at  ambient  and  low 
temperatures  as  illustrated  by  a  model  shown  in  Figure 
12.  The  point  of  bifurcation  and  load  carrying  capacity 
appear  significantly  reduced  as  the  temperature  is  de¬ 
creased.  The  postbuckling  path  also  dramatically 
changed  at  low  temperatures,  giving  the  apperance  of 
less  brittle  behavior.  In  developing  the  theory,  we 
discussed  that  this  behavior  could  be  expected  due  to 


higher  imperfection  amplitudes  induced  as  the  tempera¬ 
ture  is  decreased. 

Figure  13  shows  the  postbuckling  curves  of  the 
applied  load  ratio  K  (i.e.,  vs  the  deflection 

ratio  (8//j)  as  a  function  of  the  imperfection  ampli¬ 
tude  ratio  (8//i).  It  can  be  seen  that  room  temperature 
corresponds  to  zero  imperfection  amplitude  and  the 
plate  will  remain  flat  for  P  <  P^^.  For  a  given  imper¬ 
fection  amplitude  (due  to  thermal  load),  the  out-of¬ 
plane  deflection  occurs  even  at  no  mechanical  load.  The 
application  of  mechanical  load  will  further  increase  the 
deflection.  As  the  imperfection  amplitude  increases,  the 
above  trends  become  more  pronounced.  Therefore, 
with  reduction  of  temperature,  imperfection  amplitudes 
will  increase,  resulting  in  the  characteristic  buckling 
and  postbuckling  path  shown  in  Figure  12. 

Figure  14  shows  the  buckling  and  postbuckling 
paths  of  three  plates  tested  at-1 9,-14,  and-9.6°C.  These 
results  were  obtained  from  the  experimental  load  ratio 
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Figure  12.  Influence  of  low  temperature  on  0°  graphite! 
epo.xy  composite  buckling. 
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Figure  13.  Predicted  behavior  of  unidirectional  plates  under  buckling  load  with  initial 
imperfection  developed  from  thermal  strain. 
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and  the  corresponding  strain  data.  These  paths  lie  be¬ 
tween  the  paths  corresponding  to  imperfection  ampli¬ 
tude  ratio  (6//i)  =  0.01  and  0.3  plotted  by  using  eq  19.  If 
Mallick’s(1988)dataforT-300fiberthermal  expansion 
coefficient  of  -0.1  x  10-®  m/m_per  °C  are  used,  the 
initial  geometric  imperfection  (5//i)  from  eq  9  varies 
from  0.29  at  - 1 0°C  to  0.33  at  -20°C.  Thus,  it  can  be  seen 
that  the  theory  slightly  overestimates  the  buckling,  but 
the  trend  of  the  buckling  curves  clearly  agrees  with  that 
predicted  by  the  postbuckling  theory.  The  theory  estab¬ 
lished  that  with  lower  temperature  the  imperfection 
amplitude  is  higher. 

CONCLUSION 

The  current  work  on  the  unidirectional  graphite/ 
epoxy  laminate  is  particularly  interesting  since  such 
plates,  because  of  their  negative  thermal  expansion 
coefficient,  expand  in  the  fiber  direction  on  cooling,  but 
bend  when  constrained  without  any  mechanical  load.  A 
theoretical  treatment  of  the  problem  has  been  developed 
to  show  that  the  out-of-plane  bending  could  be  treated 
as  an  initial  geometric  imperfection.  The  postbuckling 
theory  based  on  this  initial  imperfection  model  was 
validated  by  a  series  of  experimenal  data  obtained  from 
testing  the  laminates  at  various  low  temperatures.  The 
experimental  postbuckling  curves  under  combined 
thermal  and  mechanical  loads  show  the  effect  of  in¬ 
creasing  severity  of  initial  geometric  imperfection  with 
decreasing  temperatures. 

The  implication  of  a  plate’s  bending  due  to  negative 
thermal  expansion  coefficient  is  very  significant  in 
design  considerations.  Obviously,  when  the  ambient 
temperature  drops,  such  constrained  plates  are  actually 
buckled  already  (the  structure  is  weakened  significantly) 
even  before  the  application  of  the  mechanical  com¬ 
pressive  load.  The  induced  bending  effect  at  low  tem¬ 
peratures  turns  out  to  dominate  the  stiffening  effect  in 
the  buckling  experiments.  Thus,  in  the  design  of  plates 
against  buckling  in  cold  regions,  it  is  important  to  take 
into  account  the  in-plane  constraint  causing  the  prema¬ 
ture  bending  of  the  plate  due  to  thermal  cooling  alone. 
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